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v The  capability  for  solving  the  following  two  problems,  under 
conditions  of  speed-of-sound  distribution^  a^d  other  factors  in  the  sea, 

o?  . V ^ 

is  of  well  established  importance; 

Problem  (a):  Determination  of  the  distribution  of  sound  intensity,  through- 
out the  ambient  and  extended  ocean  medium,  caused  by  a unit  point- 
source  at  a given  location. 

Problem  (b):  Determination  of  the  distribution  of  listening  sensitivity, 
throughout  the  ambient  and  extended  ocean  medium,  afforded  by  a 
receiver  at  a given  location. 


A major  component  in  the  solution  of  these  problems  by  ray  tracing  is  the 
determination  of  the  wavefront  spreading  factor  along  the  rays.  The 
present  technical  note  establishes  the  reciprocity  transformation,  for  this 
factor,  between  problems  (a)  and  (b) . 


The  Specific  Wavefront  Length 


We  focus  our  discussion  in  sound  intensity  mapping  in  a vertical 
plane.  For  a point  source  each  selected  ray  is  associated  with  an  angle, 
y , at  which  the  ray  emanates  from  the  source;  by  convention  this  angle 

s 

is  measured  clockwise  from  the  horizontal.  The  labelling  of  rays  by  their 
source-emission  angles,  Fg  , defines  Tg  wherever  the  medium  is  traversed 
by  the  rays.  It  should  be  noted  that  Tg  is  not  necessarily  single  valued 
throughout  the  medium;  in  regions,  two  or  more  different  rays  may  pass 

I 

through  every  point;  however  each  ray  passing  through  a point  belongs  to 
a particular  continuum  Of  Tg  (i.e.  family  of  rays). 

For  a family  of  rays,  in  the  vertical  plane,  the  specific  wavefront 
length,  denoted  by  L,  is  defined  at  each  point  by 
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f ' ”s  (1) 

where  n is  a unit  vector,  clockwise  normal  to  the  ray  at  that  point. 

Negative  L thus  implies  yg  is  decreasing  in  the  direction  in;  sign 
changes  in  L are  important  as  indications  of  foldings  of  the  wavefront 
as  at  caustics  or  by  reflections. 

The  differential  equation  for  integrating  L along  a ray  is  developed 
in  Reference  [l  ] as 

a2L  _ l 
as2  " ' c 

where  s is  ray  path,  C is  sound  speed  and  n is  ray  normal.  The  initial 
conditions  at  the  source  point  are 

Lg  = 0 , (dL/ds)g  = 1 (3) 

The  parameter  L may  be  integrated  along  each  ray  beginning  at  the 
source.  The  advantages  of  such  direct  integration  include  (1)  the  pre- 
clusion of  the  problem  of  measuring  specific  spreading  by  geometric 
measurement  over  finite  distances  between  diverse  rays,  (2)  the  indication, 
afforded  by  sign  changes  in  L,  of  foldings  of  the  wave  front  caused  by 
caustics  or  reflections,  and  (3)  indication,  afforded  by  the  spatial 
continuity  in  L,  of  the  adequacy  of  the  intensity  resolution  afforded  by 
the  selected  family  of  rays. 

Ill  Reciprocity 

A ray  emitted  by  a source  at  point  0 and  traced  to  a point  N is 
exactly  reversed  for  a source  at  point  N.  Ray  traces  are  reciprocal. 


52C  _idC.iL 

.2  L d s ds 
dn 
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Reference  [2]  establishes  that  the  parameter  L is  not  reciprocal 

along  a ray  trace.  Denoting  the  spreading  factor  from  a source  at  point  0 

to  the  point  N by  „ , and  for  a source  at  point  N along  the  reverse  ray 
0 , N 

to  point  0 by  LN  Q , then,  in  general, 


However  reference  [2]  suggests  that  the  parameter 


is  reciprocal.  Reciprocity  is  shown  under  conditions  where  the  sound  speed, 
C,  is  a function  of  depth  only,  and  for  other  conditions.  We  shall  show 
that  this  reciprocity  is  general,  whereby 


N.O 

C„ 


IV  Transforming  the  Governing  Equation 

We  replace  the* dependent  parameter  L by  the  dependent  parameter 
Q,  by  substituting  QC  for  L,  and  transform  Eq.  (2)  as  follows: 

csf9+23S2S  + 0sfc  . 

C ,2  1 3s  »»  W , 2 

3s  3s 


1_C  _ ^ 3C  ( 1C  IQ 

W 2 Q C 3s  5s  3s' 
3 n 


The  second  derivative  of  C along  the  curve  of  the  ray  trace  is 


developed  by  vector  calculus: 


tt  • v (tt  • V C) 


= (tt  • v)2  C + (tt  • vtt)  • v C 


(8) 


where  tt  is  the  local  ray-tangent  unit  vector,  along  increasing  ray  path,  s. 
The  first  term  on  the  right  hand  side  of  Eg.  (8)  may  be  written 


(tt  *.  V)2  C 


(9) 


where  s^  is  the  linear  measure  along  the  local  ray-tangent  axis . The 
second  term  on  the  RH  S of  Eq.  (8)  contains  the  path  curvature  term  which 
is  effected  by  refraction: 


tt  • V tt 


jn  9C 
C 3n 


(10) 


where  in  is  clockwise  normal  to  tt  and  n is  the  linear  measure  along  this 
normal  axis . Thus  Eq.  (8)  may  be  written 


32C  1_  /SC. 2 

2 C 3n 

9s* 


(ID 


Substitution  of  Eq.  (11)  into  Eq.  (7)  and  simplifications  lead  to 


I 3K  39 
Q ds  3s 


(12) 
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where 


At  this  point  we  note  a remarkable  simplification  in  form:  because  n and 
s*  are  local  orthogonal  axes,  the  terms  in  parentheses  combine  to  form  the 
Laplacian  operator  in  the  vertical  plane  of  the  ray.  Thus, 


Q V “ K + 
v 


1 « 39 

Q 9 s 3 s 

_ 


The  subscript  v denotes  that  the  Laplacian  operator  is  two-dimensional  (the 
circular  derivative)  in  the  vertical  plane  of  the  ray. 

The  initial  values  at  a source,  point  0,  where  s = 0,  are 


Q0  = 0 ; (3Q/3s)0  = l/CQ 


For  further  edification  we  replace  the  independent  ray-path  parameter 
s by  the  wave  travel  time , t , related  by 


■ s 


6s  = C 6t 


Equations  (14)  and  (15)  transform  into 


£2  = 
a«2 


- C2  V 2 K 
v 


r 


0Q  = o ; OQ/M)0  = 1 


(18) 


V Proof  of  Reciprocity 


We  write  Eq.  (17)  in  the  form 


afs 

at2 


- a Q 


(19) 


r- 


where 


C2  V 2 K 
v 


(20) 


In  the  present  section  we  prove  that  the  parameter  Q is  reciprocal: 
Beginning  with  a source  point  at  either  of  the  two  specific  end  points  of  o 
ray  path,  that  is,  with 


Qq  = 0 and  (dQ/dt)Q  = 1 , 


(21) 


the  same  value  of  Q results  in  integration  of  Eq,  (19)  to  the  other  end  of 
the  path. 

Equation  (19)  allows  a simple  and  direct  finite-difference  form; 


Q ,+Q  ,-2Q  = - a fit  Q 

vn+l  n-1  n n n 


(22) 


where  n-1,  n , and  n+1  are  subscripts  referring  values  to  three  successive 
positions  of  a wavefront  along  the  ray  path,  at  fixed  time-increments,  fit. 
Equation  (22)  may  be  written  as  a simple  recursion  formula: 
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(23) 


n+ 1 


A Q 
n n 


- Q 


n-1 


where 


A £ 2 - a 6 1 

n n 


According  to  Eq.  (21)  the  initial  values  are 


Qq  " 0 


Qj  = 6t 


(24) 


(25) 

(26) 


We  set 


Q 


1 


= 1 


(27) 


by  normalizing  the  amplitude  factor.  It  follows  from  Eq.  (23)  that 


°2  ' A1 


Q3  - Vj-1 


Q4  A3A2A1"A3'A1 


QS  ‘ A4A3A2A1-A4A3-A4A1-A2A1+1 


and  so  on  with  the  numerical  integration.  Reciprocity  is  demonstrated  for 
the  path  n = 0 to  n = 5,  for  example,  by  substituting  subscripts  in  reverse 
order  as  would  result  for  a wavefront  beginning  at  point  5 as  source  point 
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and  integrated  to  point  0.  It  can  readily  be  verified  that  Eq.  (25)  is  a 
necessary  condition  for  reciprocity. 


The  general  proof  lies  in  showing  that  the  general  solution  for  a 
wavefront  emanating  from  a source  point  0 to  a point  N , which  we  shall 
denote  by  ^ is  symmetric  in  the  position  subscripts  on  the  A’s.  That 
is,  the  solution  is  unchanged  by  replacing 

n by  N-m; 


that 


At  this  point  (and  before  turning  to  the  next  page)  those  readers 
who  are  stimulated  by  mathematical  problems  are  invited  to  pause  and 
consider  formulating  the  general  term  ^ , for  completion  of  the 
reciprocity  proof. 


Although  the  solution  is  simple  in  appropriate  formalism  the  present 
writer  admits  that  it  was  not  immediately  apparent.  The  general  solution  is 
expressed  by  the  determinant  which  has  the  elements  A^_^  / Ajj_2  ' * ’ * ’ 
, in  the  main  diagonal,  all  ones  in  the  two  adjoining  diagonals,  and  all 
other  elements  zero: 


QN  = 


Vi  i 

1 AN-2  1 


1 N-3  1 


1 A 
1 n 


(all  zeroes) 


(all  zeroes) 


1 a2  1 

1 A, 


(28) 


That  this  solution  is  symmetric  in  subscripts  is  immediately  apparent.  That 
this  is  the  solution  is  shown  by  expansion  of  the  determinant  in  terms  of  the 
minors  of  the  first  column,  which  yields 


N- 1 N-l 


- Q 


N-2 


(29) 


the  second  term  on  the  right  resulting  from  further  expansion  as  the  only 
minor  of  the  first  row. 

Thus,  in  general. 


(30) 


* 
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and 


(31) 


VI  The  Specific  Wavefront  Area  and  Reciprocity  Transformation 

A wavefront  is  a surface.  Wavefront  spreading  may  be  expressed 
in  terms  of  two  orthogonal  dimensions.  The  specific  wavefront  area  may  be 
denoted  by  E: 


E = L H 


(32) 


where  L has  been  defined  as  the  specific  wavefront  length  per  unit  radian 
of  emission  angle  at  the  source,  measured  normally  to  the  ray,  in  the 
vertical  plane  of  the  ray,  that  is,  along  m,  and  H is  the  specific  wave- 
front  length  per  unit  radian  of  emission  angle  at  the  source,  measured 
normally  to  the  ray,  in  the  horizontal,  that  is,  along  m.  The  relevant 
orthogonal  triple  of  unit  vectors  is 


m = tt  x n 


(33) 


In  reference  [2]  the  governing  equations  for  the  spreading  of  L and  H 
along  a ray  are  given  as 


I 3C  3L 
L d s 3s 


(34) 


1 iC  3H 
H 9 s 9 s 


(35) 
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where  h is  linear  coordinate  along  the  horizontal  ray-lateral  axis  defined 
by  unit  vector  cn.  These  equations  encompass  the  general  treatment  which 
permits  refraction  in  the  horizontal  as  well  as  the  vertical,  and  the  ray 
traces  are  not  confined  to  planes. 

In  preceding  sections  we  have  replaced  the  dependent  parameter  L 

by  Q: 


Q s l/C 


(36) 


We  have  transformed  the  governing  equation  into  the  form 


(37) 


where  s*  is  a linear  measure  along  the  local  ray  tangent  axis,  and  we  have 
proven  that  Q is  reciprocal,  as  defined.  In  analogous  fashion  we  replace 
the  dependent  parameter  H by  P : 

P = H/C  (38) 


The  governing  equation  transforms  into 


and  the  proof  that  P is  reciprocal  is  similar  to  that  for  Q. 


The  reciprocal  parameter  for  the  specific  wavefront  area,  E , may 
be  denoted  by  W : 


W = 


(40) 
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The  reciprocity  transformation  for  specific  wavefront  area  is 


EN,0 


(41) 
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